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Abstract 

The paper aims to study the connection between symmetries and conservation laws for the 
2D Ricci flow model. The procedure starts by obtaining a set of multipliers which generates con- 
servation laws. Then, using a general relation which connects symmetries and conservation laws 
for whatever dynamical system, one determines symmetries related to a chosen multiplier. On 
this basis, new similarity solutions of the model, not yet discussed in literature, are highlighted. 
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1 Introduction 

The concepts of symmetry, invariants and conservation laws are fundamental in the study of dy- 
namical systems, providing a clear connection between equations of motion and their solutions. 
There are many reasons for computing the symmetries and the conservation laws corresponding 
to systems of which the evolution is described by differential equations. In the recent years, a 
remarkable number of mathematical models occurring in various research domains have been stud- 
ied from the point of view of the symmetry groups theory [TH3]. The most suitable technique for 
finding classes of analytical solutions, the so-called Lie symmetry method, investigates integrability 
starting from the invariance of evolutionary equations under some linear transformations of the 
variables which define the so-called Lie group of symmetries. Conservation laws proceed from the 
conservation of physical quantities: linear momentum, mass, energy, electric charge, etc. Further- 
more, conservation laws are applied in the study of partial differential equations, as for example 
in: (i) test of complete integrability and application of Inverse Scattering Transform, (ii) study of 
quantitative and qualitative properties for pdes (Hamiltonian structure, recursion operators) and 
(iii) development of numerical methods such as finite element methods |4|5| . 

There are many interesting results on the correspondence between symmetries and conservation 
laws. As examples, an identity [6] which does not depend on the use of a Lagrangian and provides 
a relationship between symmetries and conservation laws for self-adjoint differential equations is 
derived; a direct link between the components of a conserved vector of an arbitrary partial differen- 
tial equation and the Lie-Backlund symmetry generator of the equation, which is associated with 
the conserved vector's components is obtained [7J. 

The aim of this paper is to tackle the connection between symmetries and conservation laws 
using the results mentioned above and the property of the Euler differential operator to annihilate 
whatever expression having the mathematical form of a divergence. By applying the Euler operator 
on a combination of evolutionary equations, one can determine a set of multipliers generating 
conservation laws. We will apply this procedure for the 2D Ricci flow equation, a very interesting 
model coming from the gravity theory. The Lie symmetry problem and the invariant solutions of 
this model have been already discussed [8]. Using the procedure mentioned above, new symmetries 
and solutions will be highlighted. 
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The paper is organized as follows: after this introductive section, the problem of constructing a 
conserved vector will be analyzed in the second section, by means of the conservation law multipliers' 
method and also by tacking into account the relations which impose an intimate connection between 
symmetry generators and conserved vectors. In the third section, these methods will be applied 
to the 2D Ricci flow model. A set of conservation laws corresponding to a multiplier A(x, t, y, U) 
that does not depend on the derivatives of the dependent variable U, the associated Lie symmetry 
operators and some new invariant solutions, will be obtained. Some concluding remarks will end 
the paper. 

2 Connection between symmetries and conservation laws 

For each partial differential equation (pde) or for each pde system there is a local group of trans- 
formations (called the symmetry group) that acts on the space of its independent and dependent 
variables, with the property that it maps the set of all analytical solutions to itself, and so it leaves 
the form of the equation (or system) unchanged. The widely applicable method to find the symme- 
try group associated with a pde (or pde system) is called the classical Lie method. Consequently, 
the knowledge of Lie point symmetries allows us to construct the group- invariant solutions. Two 
solutions should be equivalent if there whould be a symmetry transformation that could transform 
the one into the other. Moreover, new solutions might be obtained from the known ones: by 
applying the symmetry group to a known solution, a family of new solutions would be created. 
Let us consider a n-th order partial differential system: 

A u [x,u(x),u^{x)} = 0, v = T^j (1) 

where x = {x l , i = l,p} C RP represent independent variables, while u = {u a ,a = l,q} C R q 
dependent ones. The notation designates the set of partial derivatives of u with respect to x, 
up to the n-th order. 

Let us consider the infinitesimal symmetry operator with the general form: 

*=E^)^+X>m^ ( 2 ) 

i=l a=l 

In order to explicit the above mentioned property of classical symmetries, the Lie method [9] 
applied the criterion of infinitesimal invariance of the system ([1]) to the action of the operator ([2]). 
More precisely, the following condition ought to be imposed: 

I (b) (A v )|a=o =0 (3) 

where represents the extension of n— th order of the Lie symmetry generator ([2|). 

A conservation law for partial differential equations is a divergence expression which vanishes for 
the solutions of the pde system. The origin of conservation laws comes from physical principles such 
as the conservation of mass, momentum and energy. Although there is a well-known systematic 
method [TU] able to find conservation laws for variational pdes, the applicability of Noether's method 
is limited by the fact that many interesting pde systems are not variational. A systematic procedure 
able to find conservation laws, called the direct method, has been developed |ll|12j . The direct 
method consists of two main steps : 

(i) to determine a set of conservation law multipliers so that a linear combination of the pdes 
with the conservation law multipliers should create a divergence expression. For our given differ- 
ential system ([I]) one could define a set of conservation law multipliers {A u : M x S*( n ) — > R}, 
v = l,q, if smooth functions {Pi : M x S*( n ) — > R}, i = l,p should exist such that everywhere on 
M x SM we could have: 

A l/ [x,U(x),U {n) (x)]A u [x,U{x),U {n) (x)]= DiP l [x,U(x),U {n) {x)} (4) 
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It is known that conservation law multipliers could be found using the method |12] of the Euler 
operator. Thereby, we can solve the following system of pdes: 



K[x,U{x),U i - n) {x)}/\ u {x,U{x),U { - n) {x)] \ = 0, V p = l,q (5) 



for the unknown functions {A^, v = l,q}. 

(ii) ones having determined a set of conservation law multipliers, to find the corresponding 
fluxes in order to obtain the conservation law: 

D i P i [x,U(x),U^(x)} Mx)=u{x) = A v [x,U(x),U^(x)]A^[x,U(x),U^(x)] mx)=u{x) = 
A u [x,u(x),u ( - n \x)]A u [x,u(x),u < - n \x)] = (6) 

The key property of this set of conservation law multipliers is that their existence implies the 
existence of conservation laws. Conversely, it turns out that for non-degenerate pde systems, every 
conservation law up to equivalence must arise from a set of conservation law multipliers. Generally, 
the correspondence between sets of conservation law multipliers and equivalent conservation laws 
could be many-to-one. However, if a pde system admits a Cauchy-Kolvalevskaya form [13] and the 
sets of conservation law multipliers satisfy some mild conditions, then there is a one-to-one corre- 
spondence between each set of conservation law multipliers and each set of equivalent conservation 
laws. 

Further, a relation [7J was derived between the symmetry operator X and the components of a 
conserved vector P{. It has the form: 

X W(P i ) + PW k (e) - P k D k (C) = 0, i=T~p (7) 

where X is applied in the r— th extended form. 

The conditions (J7J) with X known, joined to the conservation law D^P 1 = 0, might be viewed 
as a system of linear partial differential equations which could be solved for the components P l , 
i = l,p of the conserved vector. Yet, ([7]) could be used as well to obtain the symmetry operators 
X associated with a given conserved vector. In the next section, the latter way will be applied. 

3 Application to the 2D Ricci flow model 
3.1 Multipliers of the model 

In this section, the relationship between symmetries and conservation laws will be applied to the 
2D Ricci flow model. Its evolutionary equation, bearing large consequences in gravity theory, has 
the form: 

u u 2 

It is a well-known equation which has been studied as a continuum limit of the Toda-type 
equation. Among the main results concerning ([8]) one may mention: (i) it could be obtained as a 
particular case of the 3D Ricci flow equation which accepts a Killing vector; (m) by linearization, it 
presents various classes [13] of solutions, these depending on the "sector" where it is defined; (Hi) 
an effective study of its Lie symmetries and invariant solutions was performed [8]. 

Now, the first task will be to obtain conservation laws for the previous model, using the multi- 
pliers' method. A multiplier A of equation ([5]) has the property : 



A 



U X y U x Uy 



u u 2 



DipifaxtyMUV)] (9) 



3 



for all functions U(t,x,y), not only for the solutions u(t,x,y) of @. 

Let us consider multipliers of the form A = A(t,x,y, U). Multipliers which depend on the first 
order and higher order partial derivatives of U could also be considered. Yet, in two dimensions, 
calculations rapidly become more complicated. Therefore, computer assisted calculations may lead 
to further conservation laws. 

The right hand side of Q is a divergence expression. Hence, this expression should vanish 
identically when applying the Euler operator which, in two dimensions, takes the form: 



Eu 



dU 



Di 



dU t 



D 



d 



+ Di 







y du„ ' ~ x du 2x 



dU x 

The determining equation for the multiplier A(t,x,y,U) is: 



+ D. 







y dU 2y 



+ D X Dy 



au. 



+ .. 



Eu 



The expansion of (jlip yields the following pde: 

'Au A 2U 



2 -^U + 



U 2 



U 



U x U y — A„r/f7 r — ArrjjU, 



yU^x 



A t + — 







(10) 



(11) 



(12) 



Because (|12j) is satisfied for all functions U(t,x,y) and because the multiplier A is chosen as not 
to depend on derivatives of U, the coefficient functions of various derivatives of U should vanish. 
Thereby, the determining system for A(t, x, y, U) is generated: 



A v = A t = A xy 







By solving the system f j 1 3 1) . one obtains the general multiplier: 

A=f(x)+g(y), V/(x), Vg(y) 



(13) 



(14) 



3.2 Conservation laws and symmetries 

The first aim of this section is to determine the conserved vector associated with the multiplier 
(|14|) . From @ and (|14|) and by performing elementary manipulations, one obtains: 



tt Uxy | UxUy 



D t [(f(x) + g(y))U}+D x 



-g(y)jf 



DiP 1 
(15) 

Consequently, when U(t,x,y) is a solution of ([5]), the previous expression leads to the following set 
of conservation laws: 



D t [(f(x)+g{y))u]+D x -g(y)^ 



■fix 



u, 



0, V/(x), Vg(y) 



(16) 



As a result, one obtains that whatever conserved vector of the Ricci flow model ([8]) with multiplier 
of the form (|14p . should have the components: 



P 1 = (/(*) + a{y))u, P 2 = -g(y)^, P 3 = -f(x)^, 

u u 



(17) 



The next aim is to illustrate how we might find the point symmetries associated with the 
conservation law (|16|) . The starting point of computations is clearly the relation ([7]) in which 
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we consider P l as known components of the conserved vector and X as an unknown symmetry 
operator. The symmetry conditions for (|17p are: 



X^iP 1 ) + P l D x {£) + P 1 D y {rf) - P 2 D X {^) - P 3 D y (v) 
x (i) (p2) + p2 Dt{ip) + P *D y (rj) - P x D t (0 - P 3 D y (Z) 
X^(P 3 ) + P 3 D t (<p) + P 3 D x (0 ~ P l D t {ri) - P 2 D x {r,) 

where the symmetry operator X is of type ([2]) and has the first extension: 









(18) 



x,y,u 



d „ d d 



dt dx 



dy 



du 



+ ■ 



dm 



+ 



: _d_ 

du. 



+ i 



,± 
du,, 



(19) 



Coefficient functions (f> x , <p y will be calculated according to the symmetry theory [9|. 

Let us suppose that A = f(x) + g(y) ^ 0. The expansion of the determining equations (fTHj) and 
separation by various monomials in derivatives of u, do generate the following differential system: 



Vx 













' , df{x) dg(y) ...... , Wc . 



dx 



dy 



U7] 



dg(y) 
dy 
df(x) 



+ [f(x) + g(y)]</> 

- g{y)[<t> - u(f> u - u(p t ] 

f(x)[(f>-u4> u -wp t ] 



(20) 



dx 

Af{*) + 9{y)]tt + <j{y)4>y 

u 2 [f(x) + g(y)]r) t + f(x)</> x 










with unknown functions </>, £, 77, (p. 

Looking to the system (|2Up . we are able to distinguish two cases for solving it. More precisely, 
the latter five equations of the system (j20|) should be solved by taking into account the following 
situations: 



Case(I) : f(x) + 0, g(y) ^ p(t), £(t, x), g(t, y),<f>(t, x, y, u) (21) 

and 

Case(II) : f(x) = 0,g(y) ^ =>• ip(t,y),^(t,x),r](t,y),(j)(t,x,y,u) or 

g(y) = 0J(x)^0^ip(t,x)^(t,x),r 1 (t,y),(t } (t,x,y,u) (22) 

3.3 Symmetries and invariant solutions for Case (I) 

Under the conditions ([2T]) . through a computational method (Maple 10 program), six solutions 
could be obtained. They correspond to the following Lie symmetry operators: 
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^1) Vai = const., V/(x) 7^ const., Vg(y) 7^ const. 



A x — ai 



dt 



(23) 



(2) V6j. = const., j = 1,3,4, Vg(y) = 62 = const. 7^ 0, V/(x) 7^ const. 

(3) Vpj. = const., j = 1, 3, 4, V/(x) = P2 = const 7^ 0, Vg(y) 7^ const. 



(24) 



. 5 a 

-PIX + » 4 — +PlM — 

ox cm 



(25) 



(4) Vsfc. = const., k = 1, 7, 7^ 2, 3, V/(x) = S2 = const 7^ 0, Vg(y) = S3 = const 7^ 



X 



in 



s 4t + + se] ^ + [( - s5 - si)y + S7] h + siu -L 



(26) 



(5) V rrii = const., i = 1,4, /(x) = (x — mi), g(y) = [— 2(y + m 2 )]( l ^ 2 ^ 



X, 



(I) 



d d d 

[m s t + m 4 ]— - m 3 [x - mi]— + 2m 3 [y + m 2 j — 

r 1 ei \ ms . . N1 c_2i d 

[m 3 t + m 4 ]—-m 3 f(x)—- — [g(y)} ( > — 



(27) 



(6) V nt = const., i = 1, 4, /(x) = (x - n 2 ), 5(2/) = (y - n 3 ) 



X 



ni 
3 

ni 
3 



t + 7T-4 

t + n 4 



5 ni(x — n 2 ) 9 n\{y — n 3 ) d 
_ 



5 



3 <9x 3 
_<9 _ ni/(x) 9 _ nig(y) d 
5t 



3 cte 3 dy 1 ~du 



dy ' niU du 
d 



(28) 



It is important to remind some results [8]: 

(i) by applying the classical symmetry approach, the general Lie operator for the Ricci flow 
model jHl) is obtained in the form: 



^ = (c 1 t + C2)| + C(x)A +r?(y) | +n 



d£{x) _ dn{y) 
dy 



dx 



d_ 

du 



(29) 



with ci, C2 arbitrary constants and £(x), n(y) arbitrary functions. 

(ii) in the linear sector of invariance (when £(x), n(y) are considered as having linear forms), 
the Ricci flow model admits a 6— parameters family of Lie operators which generates the following 
6 independent symmetry operators: 



d d Tr d Tr d Tr d 
Vi = t— + u—, V 2 = — V 3 = —, V 4 = — , 
at Ou at ox oy 

d d d d 

V $ = x a W ^~' V 6 = 2/^ U— 

ox ou oy ou 



(30) 



Remark 1 : The Lie point symmetries associated to the conservation law (|16p . which have the 
expressions (f23|) - ([28|) could be expressed as a linear combination of the independent operators ([30]) . 
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These expressions are: 

X[ I] = ai V 2 , = b 3 V 2 + b A V 4 - hV 6 , X { 3 1] = p 3 V 2 + P xV 3 - Pl V 5 (31) 

X { A I] = s 4 V 2 + s 6 V 3 + s 7 V 4 + s 5 V b -{sx + s 5 )V 6 (32) 

xjp = m 3 Vi + m 4 V 2 + mim 3 V 3 + 2m 2 m 3 Vx - m 3 V 5 + 2m 3 V G (33) 

XP = y V X + n 4 V 2 + nm 2 V 3 + n x n 3 V A - ^ V 5 - ^ V 6 (34) 

The next objective of this section is to determine the invariant solutions of the analyzed model, 
generated by the Lie operators ([23]) - ([28]) . The similarity reduction method [15] will be applied. 
Consequently, the operator x[^ from (J23J) has the characteristic equations: 

dt dy dx du 

ax~ = ~0~ = T = ¥' ( ' 

By integrating these equations, one obtains three invariants: 

/« = x, = y, if = u. (36) 

By designating the invariant 1^ = Hx(x, y) as a function of the other two and then inserting it to 
the Ricci flow equation ([8]), the similarity reduced equation takes the form: 

(Hx) x (Hx) y -Hx(Hx) xy = 0. (37) 

The solution of this equation is: 

Hx(x,y) = hx(x)h 2 (y), V hx(x), Vh 2 (y) (38) 
Thereby, the invariant solution corresponding to the operator X± has the final form: 

4 /} (x,y) = hx(x)h 2 (y), V hx(x), Vh 2 (y) (39) 

Through a similar algorithm, the following results could be obtained: 
(1) The Lie operators (|24[) . (|25p . (|26|) do generate, respectively, the following similarity variables: 

hy-bi fh,\ pxx-pi fpx 
v 2 = x, w 2 = exp [ —t ; v 3 = y, w 3 = exp — t 



h \h J " Px \P3 

s$x + sq ( s 5 \ {s 5 + sx)y-s 7 fs 5 + sx,\ ,. n . 

vx = exp 1 , w 5 = ■ exp 1 (40) 

s 5 V s 4 / s 5 + sx \ s 4 J 

and, through similarity reduction, evolutionary equations will look like: 

bxw 2 (H 2 ) 2 (H 2 ) W2 + bx {H 2 f - 
- b 3 H 2 (H 2 ) V2W2 + b 3 (H 2 ) V2 (H 2 ) W2 = (41) 

2 try \ 1 „ I tj \3 



Pl w 3 (H 3 Y(H 3 ) m +px (H, 



:-;) 



-p 3 H 3 (H 3 ) V3W3 +p 3 (H 3 ) V3 (H 3 ) m =0 (42) 
■■ 5 vx {H A f (Hx) V4 - (s 5 + Sl )wx (Hx) 2 (Hx) W4 - 

Sl (Hx) 3 + s 4 Hx(Hx) V4Wi - sx {H 4 ) V4 (Hx) W4 = (43) 
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(2) The invariant solutions of the model, associated to X 2 \ X§ respectively xjjP , become 
stationary: 

4 1} (^y) = j-^-Fzix), 4 I \x,y) = ^± G 3 (y), VF 2 (x), V G 3 (y) (44) 

W f « 5 * + ^ V (* 5 + s x )y - s 7 \ (*p-0/("+i) 
«4 O^y) = ft ; , Vgi,g 2 = const. (45) 

(3) By integrating characteristic equations, the Lie operators (|27p and (|28p do lead to the following 
new variables: 

( u y + m 2 

v 5 = m 4 x + m 3 (x - mi i, w 5 = — ^ 

[m 3 t + m 4 y 

vq = x{n\t + 3114) — nin 2 i, wq = y{n\t + 3114) — n\n 3 t (46) 
Through the similarity reduction method, they generate the reduced differential equations: 
m 3 (v 5 - mim 4 ) (H 5 ) 2 (H 5 ) Vs - 2m 3 w 5 (H 5 ) 2 (H 5 ) Ws - 

- H 5 (H 5 ) V5W5 + (H 5 ) V5 (H 5 ) W5 = (47) 
ni(«6 ~ 3n 2"- 4 ) (H 6 f (H 6 ) V6 + m(w 6 - 3n 3 n 4 ) (H G ) 2 (H 6 ) We + 

+ 3n x (H 6 f - H 6 (H 6 ) V6W6 + (H 6 ) V6 (H 6 ) W6 = (48) 

(4) The invariant solutions of the model, corresponding to z] , respectively to Xq are nonsta- 
tionary. They have the expressions: 

uip(t,x,y) = Pi (mim4 — v^)w^ 2 ^ , V /t> 1} p 2 = const. 

u£\t,x,y) = 7 1 (u 6 - 3n 2 n 4 ) 72 (-u> 6 + 3n 3 n 4 )~ (3+72) , V7 1 ,7 2 = const. (49) 
where invariants V5, 105, t>6, u>6, expressed in the terms of original variables, are provided by (|46|) , 

3.4 Symmetries and invariant solutions for Case (II) 

By solving the latter five equations of the determining system (|20p under the conditions ()22|) . the 
following Lie symmetry operators will be generated : 

^> = 4 + - (^«) I;, c 2 = con,,, V«.) (50) 



X 2 = c 2 -+ V (y)--^— u)-,c 2 = const. V V (y) (51) 



Remark 2: The forms (|50p , ()5ip represent particular cases of the general operator (|29p . They could 
be obtained by imposing c\ = 0, n(y) = 0, respectively ci = 0, £(a?) = in (f2§j) . 

Let us obtain invariant solutions for the 2D Ricci model, associated to operator (|50p . Similar 
results will be derived for operator ([5T]) . 

The function -u = *$>(t,x,y) is an invariant solution generated by (|5Q|) . providing that: 

x{ /7) [u - *(t,x,j/)] | u= * = (52) 
The previous condition has the equivalent form: 

+ ca^t + e(x)* a = (53) 
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The general solution of f[53|) is: 



F (v w) f Co 

*(t, x, y) = l\> , v = y, w = t- TT^ dx ( 54 ) 

with arbitrary functions £(cc) and F (v, w) . 

The substitution of (|54f) into the 2D equation ([8]) yields a reduced differential equation of the 
form: 

F 2 F W + c 2 FF vw - c 2 F v F w = (55) 

which admits as a solution an arbitrary function F±(v). 

Consequently, by taking into account (|54p . the invariant solution associated to operator (|50p 
becomes: 

u[ II) (x,y) = ^, Viq(y), V£(x) (56) 
Due to similar reasons, the invariant solution attached to operator (|5ip will be: 



u ^\ x ,y) = ^l, VF 2 (x), Vr?(y) (57) 

my) 



4 Concluding remarks 

In literature the problem of constructing a conserved vector is generally approached by means of the 
direct method of solving DiP 1 = 0, for a given equation with no recourse to symmetry properties. 
It usually involves ad hoc assumptions able to simplify the procedure. For a known symmetry 
operator X, determining the conserved vector P = (P l ), i = l,p whould become simpler if the 
relation ([7]) which connects symmetries and conservation laws, should be added to the conservation 
law DiP 1 = 0. Conversely, the special relation ([7]) could also be used to obtain the symmetry 
generators X associated with a given conserved vector P. The latter approach is applied in this 
paper for the 2D Ricci flow model. 

The central objects, in the study of conservation laws, are the sets of conservation law mul- 
tipliers. The key property of these multipliers is that their existence implies the existence of 
conservation laws. The main results of this paper consist in obtaining: (z) the set of multipliers 
([11]) which depend up on two arbitrary functions f(x) and g(y); (ii) the set of local conservation 
laws (fT6|) associated to the mentioned multipliers; (in) two sets of symmetry operators (f23|) - ([28j) 
and (|50p , (|5ip , by making use of the the tight relationship ([7]) between symmetries and conservation 
laws and by solving the determining differential system (|20p in two different cases. In the first case, 
each of the Lie symmetry generators could be expressed as a linear combination of the independent 
operators (|30p . obtained [8J for the linear sector of invariance of the Ricci flow model. In the sec- 
ond case, both operators (f5Pj) and (|5Tj) . could be identified as particular forms of the general Lie 
operator [8] which depends on two arbitrary functions; (iv) six stationary and two nonstationary 
invariant solutions of the Ricci model, associated to the mentioned symmetry operators are com- 
puted. Among, them, some new similarity solutions, not yet given in literature, were highlighted. 
They are listed in (g5J, (06]), (09|), ([56]), floTD . 

These results have been obtained by assuming a conservation law multiplier which does not 
depend on the derivatives of the variable U(t,x,y). The multipliers which depend on the first 
order and higher order partial derivatives of U may lead to further conservation laws. They will 
be studied in forthcoming works. 
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